Draft version March 20, 2009 

Preprint typeset using I^T^X style cmulatcapj v. 10/09/06 



DISSECTING THE RED SEQUENCE— II. STAR FORMATION HISTORIES OF EARLY-TYPE GALAXIES 

THROUGHOUT THE FUNDAMENTAL PLANE 

Genevieve J. Graves 1 , S. M. Faber 1 , & Ricardo P. Schiavon 2 
Draft version March 20, 2009 

ABSTRACT 

This analysis uses spectra of ~16,000 nearby SDSS quiescent galaxies to track variations in galaxy 
star formation histories along and perpendicular to the Fundamental Plane (FP). We sort galaxies 
by their FP properties (a, R e , and I e ) and construct high S/N mean galaxy spectra that span the 
breadth and thickness of the FP. From these spectra, we determine mean luminosity-weighted ages, 
[Fe/H], [Mg/HI, and [Mg/ Fe] based on single stellar population models using the method described in 



Graves & Schiavon (2008). In agreement with previous work, the star formation histories of early-type 
galaxies are found to form a two-parameter family. The major trend is that mean age, [Fe/H], [Mg/H], 
and [Mg/Fe] all increase with a. However, no stellar population property shows any dependence on 
R e at fixed a, suggesting that a and not dynamical mass (Md yn a cr 2 R e ) is the better predictor 
of past star formation history. In addition to the main trend with a, galaxies also show a range of 
population properties at fixed a that are strongly correlated with surface brightness residuals from 
the FP (Alog/ e ), such that higher surface brightness galaxies have younger mean ages, higher [Fe/H], 
higher [Mg/H], and lower [Mg/Fe] than lower-surface brightness galaxies. These latter trends are a 
major new constraint on star-formation histories. 

Subject headings: galaxies: abundances, galaxies: elliptical and lenticular 



1. INTRODUCTION 

Early-type galaxies are known to comprise a two- 
parameter family in galaxy structure. In the three- 
dimensional (3-D) parameter space of central stellar ve- 
locity dispersion (a), galaxy effective radius (R e ), and 
mean central surface brightness (I e ), early type galaxies 
occup y a 2-D plane, known as the Fundamental Plane 
(FP; |Faber et al.[|1987| [Dressier et"aL]|1987| |Djorgovski| 
|fc Davis|1987p . Th e FP is relatively cons tant throughout 
the local universe (|j0rgensen et a l.|1996|). There may b e 

Hernardi 



et al.|2003a| , 

nts being ofi- 



some variation with environment ( I 
consistent with the FP in dense environments 
set toward lower I e than the FP in the field. 

There is clear evolution in the FP with redshift, 
which is often interpreted as evolution in the stellar 
mass-to-light ratio (M/L). For field galaxies, most au- 
thors find consistent results, with M/L ev olution of 



d\o%(M/L B )/d z 0.7 (Treu et al.||2005bl and refer- 
ences therein). Treu et al.| Q2005b"| also hnd increased 
scatter in the FP at higher redshift, with most evolu- 
tion seen in the low-mass end; while the high-mass FP 
appears nearly unchanged since z ~ 1 other than ex- 
pected passive evolution in M/L, low- mass galaxies are 
still "settling" onto the FP since z ~ 1 (pan der Wei et al. 



2004 Treu 



et al.||2005a|b| F 
il FP should Hi 



van der Wei et al. 



The local FP should then be expected to include at 
least some galaxies that are relati vely recent a r rivals 



Thi s is supported by th e resul ts of jForbes et al. (19981 



and jTerlevich fc Forbes ( |2002[ ) , who hnd that residuals 



from the FP correlate with age such that galaxies offset 
to higher (lower) surface brightness have younger (older) 
ages than those occupying the FP. The thickness of the 
FP may therefore be an age sequence, with the most re- 
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cent arrivals lying at higher surface brightness than the 
midplanc. This thickness may reveal a (narrow) third di- 
mension to the structural parameters of early-type galax- 
ies. 

In contrast to the 2-D structural space occupied by 
early-type galaxies, most studies of their stellar pop- 
ulations have stressed only a 1-D space. Many pre- 
vious authors have presented correlations between de- 
tailed stellar population properties and various measures 



of galaxy "size", including a (e.g., jTrager et al.|[2000| 
Thomas et al.|2005|[Ne"lanet al.|20 05; Smi thet al.|2007a| , 
stellar mass M* (e.g., Gallazzi et al. 2005), and lumi- 
nosity (L) (e.g. , Kuntschner & Davies 1998; Terlevich 
& Forbes 2002 1 . However, considerable evidence sug- 



gests that early-type galaxy star formation histories in 
fact c omprise a two parameter family. |Worthey et al.| 
( |1995[) demonstrated an anti-correlation between age and 
metallicity among galaxies with similar a and showed 
that it helped explain the narrown ess of the color-cr and 



FP relations. Trager et al. (|2000[) quantified this anti- 



correlation, which they termed the "metallicity hyper- 
plane." The anti-correlation has been corroborated re- 



cently by |Smith et al. (|2008bJ. |Thomas et al.| ( |2005| | also 
found a range of ages for galaxies at fixed a although with 
no associated metallicity variation. 

This is the second paper in a four-part series which 
demonstrates conclusively that early-type galaxies do 
span a 2-D space in stellar population properties, im- 
plying that their star formation histories must also oc- 
cupy a two-parameter space. The goal of this series is to 
connect stellar population properties of early-type galax- 
ies to their structural parameters in order to understand 
how galaxy structural evoluti on and mass assemb ly are 
related to star formation. In Graves et al. (2009 here- 



after Paper I) , we explored the variations in stellar pop- 
ulations along the color- magnitude relation of early-type 
galaxies and showed that galaxy stellar populations not 
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only vary systematically with a, but also vary with resid- 
uals from the a-L relation. 

In this paper, we continue that approach but map out 
stellar populations along and through the FP rather than 
the a-L relation. The result is to parameterize the 2-D 
space of stellar population variables in terms of the FP 
variables a, R e , and I e . Briefly looking foward, Paper 
III presents a clear parameterization of the 2-D family 
of early- type galaxy star formation histories and inter- 
prets the observed stellar population results as differ- 
ences in the onset epoch and duration of star formation in 
galaxies. Paper IV compares various stellar (M*/L) and 
dynamical (M^yn/L) mass-to-light ratios for the sample 
galaxies and demonstrates that the observed stellar pop- 
ulation effects are not adequate to explain the variation 
in Mdyn/L along or perpendicular to the FP. 

The paper is organized as follows. Section [2] briefly 
reviews the selection criteria used to identify the early- 
type galaxy sample for this series of papers. In S|3] we 
divide up the FP into bins of similar galaxies, stack the 
spectra to obtain high S/N mean spectra for each point 
in FP-space, and compare the spectra to models in or- 
der to determine the mean stellar population properties 
along and through the FP. Section [4] presents these re- 
sults. A major new result is that stellar populations, and 
therefore galaxy star formation histories, appear to be in- 
dependent of R e , at fixed a. Section [5] briefly discusses 
some possible explanations for why galaxy star forma- 
tion histories do not depend on R e at fixed a. Finally, 
^summarizes our conclusions. 

2. SAMPLE SELECTION 

The sample of galaxies used in this work is identical to 
that of Paper I. A brief summary of the sample selection 
is included here. The sample consists of ~16, 000 galaxies 
from the Sloan Digital Sky Survey (SPSS; [York et ah 



2000) spectroscopic Main G alaxy Survey ([Strauss eta L 
2002)) Data Release 6 (DR6; |Ad^an-McCarthy||2008[ ).' 
Galaxies are chosen from a limited redshift range (0.04 < 
z < 0.08) and are selected to be quiescent galaxies based 
on the following criteria: 

1. Galaxies have no detected emission in cither the 
Ha or [On]A3727 lines. 

2. Galaxies have centrally concentrated light profiles, 
as determined by the ratio of the Petrosian radii 
that enclose 90% and 50% of the total galaxy light 
(R90/R50 > 2.5 in the i-band). 

3. The likelihood of a de Vaucouleurs light profile fit 
is at least 1.03 times larger than the likelihoood of 
an exponential light profile fit. 

Paper I showed that, although no explicit color selection 
has been applied, these criteria define a sample of galax- 
ies that populate the red sequence in the color-magnitude 
diagram, with very few color outliers (see Figure 1 of Pa- 
per I). Remaining color outliers are excluded from this 
analysis, as described in Sj3] 

By requiring the sample galaxies to have no detectable 
emission, we remove four types of galaxies from the sam- 
ple: actively star-forming galaxies, Seyfert hosts, low ion- 
ization nuclear emission- line region (LINER) hosts, and 



so-called "transition objects" (TOs) whose spectra con- 
tain emission from both star-formation and active galac- 
tic nuclei (AGN). In galaxies with ongoing active star- 
formation (star- forming galaxies and TOs) , a small frac- 
tion of very young stars contribute disproportionately to 
the integrated light of the galaxy. These galaxies will 
therefore have significantly lower M*/L than quiescent 
galaxies of comparable mass, which will change their lo- 
cation in FP space. These galaxies constitute ab out 11% 
of morphological ly-selected early type galaxies (ISchaw- 
inski et al.||2007|). Excluding these galaxies removes the 



"youngest" early type galaxies from the sample, leaving a 
sample of galaxies whose M/L values and positions in FP 
space are not strongly biased by a small sub-population 
of very young stars. 

Removing galaxies with emission due to AGN activ- 
ity (i.e., Seyferts and LINERs) is a less obvious choice, 
as these galaxies likely lack the very young stars that 
bias M/L estimates in actively star- forming galaxies. 
However, host galaxies of strong AGN have stellar light 
profiles i ntermediate between e arly-type and late-type 
galaxies ( Kauffmann et aL]|2003[ ), as indicated by values 
of Rgo/n-50- Furthermore, there is evidence that these 
galaxies have systematically different stellar populations 
from quiescent ear ly type galaxies of the same mass, 
with both Seyferts jKauffmann et al.|2003j ) and LINERs 
( Graves et al.|[2007| showing younger stellar population 
ages. In light of these differences, we defer a FP analysis 
of these galaxies to future work. 

Defined in this way, the galaxy sample presented here 
contains bulge-dominated galaxies, some of which har- 
bor significant disk components. Figure 2 in Paper I 
shows a selection of the sample galaxies illustrating this. 
Classical bulges of disk galaxies a re known to lie on the 
FP defined by early-type g alaxies (Fisher fc Drory|2008 



Kormendy fc Fisher||2008| and are therefore relevant to 
this study. Because the SDSS images are fairly low- 
resolution, high quality automated bulge-disk decompo- 
sitions are challenging. Assessing the quality of the de- 
compositions and interpreting the output is a compli- 
cated topic and beyond the scope of the analysis pre- 
sented here. We are engaged in ongoing work to asses 
and compare automated bulge-disk compositions with 
by-eye morphologies (Cheng et al., in preparation). We 
then intend to revisit the topic of morphological varia- 
tion among quiescent galaxies and its effect on the FP 
and stellar population. 

The FP explored in this work is defined by the galaxy 
central velocity dispersion (a), galaxy effective radius 
(R e ), and galaxy effecti v e surface br ightness (7 e ). Fol- 



lowing U0rgenserietaL (1996 20061), values of R e are 



derived using r 1 / 4 de Vau couleurs fits to the galaxy light 
profiles (i.e., Sersic 1968 fits with n = 4). It is likely 
that not all sample galaxies will have true Sersic n = 4 
profiles and that this will affect the derived FP correla- 
tions. However, it is not clear that a physically meaning- 
ful comparison can be made between values of R e derived 
from fitting profiles with different shapes. A separate 
study of the FP relationships presented here as a func- 
tion of Sersic n will be included in the FP morphology- 
dependence study discussed above. 

In this work, velocity dispersions are taken from the 
SDSS DR6 catalog and corrected to a "central" |i? e 
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aperture (see Paper I for details) . Effective radii are from 
r-band de Vaucouleurs fits to the galaxy light profiles, 
converted to physical units using standard ACDM cos- 
mology with Q\ = 0.7, Om = 0.3, and ho — 0.70. Sur- 
face brightnesses are computed in the Johnson V band as 
I e = Lv/2ttRI, where Ly is the total V band luminosity 
of the galaxy from a de Vaucouleurss fit to the light pro- 
file. The K-correction from observed SDSS ugriz pho- 
tometry to rest-frame U-band is performed using the IDL 
code kcorrect v4.1.4 flBlanton et"al~1|2003| ). The de Vau- 
couleurs radii and luminosities are corrected for known 
problems with sky-subtraction around bright galaxies in 
the SDSS photometric pipeline, as described in Paper I. 
Luminosities are also corrected for Galactic extinction 
using the extinction values from the SDSS photometric 
pipeline. 

The data used to identify the sample and to compute a, 
R e , and I e for each galaxy include the foll owing parame- 
ters f rom the NYU Value-Added Catalog (Bla nton et al.| 
20051 vers ion of the SDSS DR4 (|Adelman-McCarthy| 
redshift (z), 



et al 



20061): redshift (z), de Vaucouleurss photometry 
in the ugriz bands, and Petrosian radii (i?go and i?so). 
These are supplemented with the following parameters 
from the SDSS DR6 Catalog Archive Server 3 : velocity 
dispersion as measured in the SDSS 3" fiber (a fib), r- 
band de Vaucouleurs radius (r<j e y) and axis ratio (a/b), 
and the likelihoods of de Vaucouleurs and exponential 
fits to the galaxy light profiles. The non-detection of 
emission line s is determined usi ng the emission line mea- 
surements of Yan et al. (2006) by requiring galaxies to 
have emission line fluxes below a 2a detection in both 
the Ha and the [On]A3727 lines. The spectra used in 
this analysis are downloaded from the SDSS DR6 Data 
Archive Server 4 . 

The galaxy radii are critical measurements for this 
analysis, entering the FP determinations both directly 
through R e and indirectly through the calculation of I e . 
We have attempted to correct for known problems in the 
SDSS photometry pipeline (see Paper I) but problems 
may remain. To check the validity of the de Vaucouleurss 
R e available in DR6, we compared these values to effec- 
tive radii measured indepe ndently using the photometry 



package GALFIT v.2.0.3c riPeng et al |2002[) kindly pro 
vided by A. van der Wei. The GALFIT radii were mea- 
sured from Sersic fits to the galaxy r-band light profile 
with n — 4 to match the SDSS de Vaucouleurss photom- 
etry In comparison to the GALFIT radii, the SDSS radii 
arc slightly underestimated for the largest galaxies (low 
by ~ 0.07 dex) and slightly overestimated for the small- 
est galaxies (high by ~ 0.1 dex). In addition to these 
systematic variations, there is scatter of ~ 0-05 dex be- 
tween the two estimates of logi? e . In |3.1| galaxies are 
sorted into bins based on R e , with bin widths of 0.2 dex 
in log R e . Uncertainties on the order of 0.05 dex in log R e 
should therefore have only modest effects on bin assign- 
ments. To make certain that our results are not strongly 
affected by possible errors in the DR6 values of R e , we 
have performed the entire analysis presented in this pa- 
per using both the SDSS pipeline values of R e and the 
GALFIT fits to R e and compared the differences. All 
the results are qualitatively identical. We have chosen 

3 http://cas.sdss.org/dr6/en 

4 http: //das.sdss.org/DR6-cgi-bin/DAS 
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Fig. 1. — The thickness of the FP in the I e dimension. Here, 
the measured value of log/ e for each galaxy in the sample is 
plotted against the predicted logI e in terms of log a and log R e 
from the best fit to the FP. We divide the FP into "low-SB" 
(-0.3 < Alog/e < -0.1), "midplane" (-0.1 < Alog/ e < 0.1), 
and "high-SB" (0.1 < Alog/ e < 0.3) slices based on the surface 
brightness residuals, as indicated by the solid lines. 

to present the results based on the SDSS DR6 pipeline 
photometry as these values are more easily accessible to 
the general astronomical community. 

3. DECONSTRUCTING THE FUNDAMENTAL PLANE 

The overarching goal of this project is to understand 
the relation between galaxy star formation histories 
(SFHs) and their present-day structural properties, as 
categorized by their location in FP- space. To study 
systematic variations in SFHs with structure, galaxies 
are sorted into bins based on their location on the FP, 
then average spectra are constructed from the individ- 
ual galaxies in each bin. The average spectra have very 
high signal-to-noise (S/N), allowing accurate stellar pop- 
ulation modelling of the mean SFH of galaxies at each 
point in FP-space. This section defines the binning and 
describes the stellar population modelling process. 

3.1. Binning in Fundamental Plane Space 

In addition to studying stellar population trends over 
the FP, we also want to investigate variations through the 
thickness of the FP. We thus want to divide galaxies into 
bins in a three dimensional (3-D) space. There are two 
obvious choices for binning coordinate systems. The first 
is to use two orthogonal vectors within the FP, together 
with a third vector d efined perpendicular to the FP, sim- 
ilar to the K-space of |Bender et al. ( 1992 1. The second is 
to use the observed structural parameters log a, logi? e , 
and log I e directly, despite the fact that they are not truly 
orthogonal. We have chosen the latter approach because 
the resulting trends yield more readily to interpretation. 

The parameters log a and log R e fortunately describe a 
relatively face-on view of the FP. It is furthermore attrac- 
tive to choose them as fundamental binning parameters 
because neither depends explicitly on stellar population 
properties (in contrast to I e ) and because a combination 
of a and R e gives an estimate of the dynamical mass 
(Mdy n ) of a galaxy, thus this binning scheme groups to- 
gether galaxies of similar total mass. The FP by defini- 
tion has only two independent parameters. Once log a 
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Fig. 2. — Bin definitions for the purposes of sorting and stacking galaxies across the FP. Gray points show all galaxies in the sample, 
while black lines indicate the bin definitions. The three panels represent the three slices in surface brightness residuals, defined in Figure 
^ Galaxies are further divided into 6 bins in logcr and 5 bins in logi? e within each slice in Alog/ e , as illustrated. Five of the 90 bins 
do not contain enough galaxies to produce robust mean spectra and are therefore excluded from our analysis. These are indicated by the 
shaded bins. In each of the three surface brightness slices, the galaxies occupy somewhat different parts of log c-log R e space, indicating 
some curvature in the observed FP for this sample. 



and log R e have been set, the mean value of log I e for 
that bin is fixed. 

We use the log I e dimension to explore SFH variations 
through the thickness of the FP. Treating log I e as the 
dependent variable, we fit the FP for log I e as a function 
of log a and log R e using a simple least squares fit that 
minimizes residuals in the log I e direction (implemented 
as the IDL routine sfit.pro) and find the relation 

logJ e = 1.16 logcr- 1.21 log R e + 0.55, (1) 

where a, R e , and I e are in units of km s _1 , kpc, 
and Lq pc -2 , respectively. The thickness of the FP 
can then be described as Alog/ e , where A indicates 
the difference between the observed value of log I e and 
the expected value of log/ e given the observed logcr 
and log R e . Figure [T] shows the observed log I e plot- 
ted against the best fitting FP relation given by equa- 
tion [T] By defining cuts in A log I e , the FP is sliced 
into three layers: "low-SB" (-0.3 < Alog/ e < -0.1), 
"midplane" (-0.1 < Alog/ e < +0.1), and "high-SB" 
(+0.1 < Alog/ e < +0.3), as shown in Figure [I] Here, 
the terms "low-SB" and "high-SB" are relative terms 
which describe the value of I e with respect to the typical 
value of I e for galaxies of the same cr and R e . 

More sophisticated statistical methods can be used to 
fit the FP relation, but the assignment of galaxies to 
bins is relatively insensitive to the exact equation used 
to describe the FP. The choice of log I e as the dependent 
parameter in the binning process is motivated by three 
things: the first is the previously stated fact that I e (and 
not R e or a) is expected to vary most strongly with star 
formation history. The second is that a and R e can be 
combined to estimate M^ yn for a galaxy, thus produc- 
ing groups of galaxies of similar total mass. The third is 
that Paper I demonstrated that, at fixed cr, stellar popu- 



lations vary substantially with A log Ly, which is related 
to Alog/ e . 

In each of the three surface-brightness slices, galaxies 
are divided by a 6 x 5 grid in log cr and log R e , as shown in 
Figure[5J The grid is identical for all three slices. Sorting 
galaxies in this 3-D space of log cr-log R e -A log I e results 
in 90 different bins of galaxies. Within each bin, galax- 
ies share similar values of these three parameters, and 
therefore also share similar values of Md yn , total lumi- 
nosity (Ly), and mass-to-light ratio {M,i yn /L). Within 
each bin, galaxies with colors more than 0.06 mag 5 from 
the mean g — r color of the bin are excluded so that the 
stellar population properties measured in each bin are 
not biased by a small number of interlopers. The me- 
dian values of logcr, logi? e , Alog/ e , and log/ e for each 
of the 90 bins are listed in Table 1, along with the total 
number of galaxies in each bin. 

Interestingly, Figure [2] shows that the galaxy distribu- 
tion in log cr-log R e space is somewhat different in the 
various slices, with the high-SB slice showing a tail of 
high cr, high-i? e galaxies that are not present in the 
low-SB slice. This occurs because the FP of our sam- 
ple galaxies has a small degree of curvature. This curva- 
ture is not necessarily real and may be due to selection 
effects. We have not attempted to correct for selection 
biases in this analysis, which spans a limited range in 
rcdshift and is therefore nearly volume- limited. Because 
galaxies are binned by surface brightness, selection ef- 
fects should result in low-SB bins being incomplete but 
should not substantially bias the mean spectra. 

Of the 90 bins described here, six bins do not contain 
enough galaxies to produce a reliable stacked spectrum 

5 This is the typical color spread for galaxies with the same a. 
See section §4.1 of Paper I for details. 
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(i.e., contain 5 or fewer galaxies), leaving a total of 84 
stacked spectra for the analysis. Spectra are combined 
using a method that rejects outlier pixels, masks bright 
skylines, weights all galaxies equally, and smoothes all 
spectra to the same resolution before combining (see §4.2 
of Paper I for details). The result is a high S/N mean 
spectrum for each bin on the FP, as well as a correspond- 
ing error spectrum. 

The bins are not evenly populated. As can be seen 
in Figure [2] and Table 1 , some of the bins contain many 
hundreds of galaxies, while others contain a dozen or 
so. The typical S/N of the individual stacked spectra 
therefore varies substantially between bins. It is worth 
bearing in mind that the error bars are not uniform for 
the stacked spectra. In much of the subsequent analysis, 
we will indicate data derived from the lower S/N spectra 
and present typical error bars separately for the low-S/N 
and high-S/N data. 

3.2. Lick Indices and Stellar Population Modelling 
In each of the stacked spectra, we measure the full 



set of Lick indices as defined in Worthey et al. (19941 



and Worthey & Ottaviani (1997[). These include Balmer 



lines H/3, ti'fF, and Hd^ (as well as broad versions of the 
bluer Balmcr lines, Hja and H5a), a set of Fc-dominated 
lines (Fe4383, Fe4531, Fe5015, Fe5270, Fe5335, Fe5406, 
Fe5709, and Fe5782), and lines that are also sensitive 
to element abundances other than Fe (Mgi, Mg2, Mg b, 
CN 1; CN 2 , C 2 4668, Ca4227, Ca4455, NaD, Ti0 1; and 
TiOa). Index measurements are performed by the au- 
tomated IDL code Lick JEW, which is part of the pub- 
licly availa ble 6 EZ-Ages code package (Graves & Schi- 



avon 



2008[ ). Lick_EW also computes errors in the Lick in- 
dex absorption strengths from the associated error s pec- 
tra following the formalism of |Cardiel et al. (19981. It 
should be noted that the SDSS 3" spectral fibers sample 
roughly 2/3 R e for the redshift range of our sample. In- 
dex measurements and derived abundances in this work 
are therefore more nearly comparable to literature values 
measured within R e or R e /2, rather than central or R e /8 
values. Because many early- type galaxies show strong 
radial mctallicity gradients, central index strengths and 
abundances tend to be higher than those measured from 
a substantial fraction of the galaxy light. 

All of the stacked spectra have been smoothed to 
match the highest-er galaxies in the sample (300 km s _1 ) 
in order to compare all galaxies at the same effective 
spectral resolution. The combined smoothing due to a 
and the SDSS native resolution is lower than the resolu- 
tion at which the Lick indices are defined (~ 210 km s _1 
at SDSS resolution); thus the linestrength measurements 
must be corrected to bring them onto the Lick system. 
We do th i s usin g corrections computed in Table A2a of 
Schiavon (20071. The stellar population models used to 
interpret the line strength measurements are based on 
flux calibrated spectra, and the SDSS spectra are them- 
selves flux-calibrated; thus we do not apply a zeropoint 
shift to bring our index measurements fully onto the Lick 
system. Because the data and the models are both flux- 
calibrated, zeropoint offsets should be negligible. 

To convert the set of line strength measurements into 
stellar population mean ages and abundances, the Lick 

6 http:/ /www. ucolick.org/~gravcs/EZ_Agcs. html 
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Fig. 3. — Lick index measurements for the stacked spectra, com- 
pared with a grid of stellar population models. The stellar popula- 
tion models shown here are computed using the average abundance 
pattern for the sample galaxies ([Mg/Fe] = +0.19, [C/Fe] = +0.17, 
[N/Fe] = +0.05, [Ca/Fe] = +0.01, with [Na/Fe] = [Si/Fe] = [Ti/Fc] 
= [Mg/Fe], and [Cr/Fe] = 0.0. See |Graves fe Schiavon|2008| for de- 
tails). Solid lines connect models at constant [Ke/HJ, while dotted 
lines connect models at constant age. A typical statistical error for 
the index measurements is shown in the top left. Lines of constant 
age (metallicity) do not run exactly horizontal (vertical) because 
the H/3 and (Fe) indices are sensitive to both age and metallicity 
(the "age-metallicity deneracy"). As is clear from the model grids, 
this degeneracy can be broken using a combination of Lick indices. 
The non-orthogonality of the model grids results in correlated er- 
rors in the age and [Fe/H] determinations for the stacked spectra, 
as discussed in 94.31 



indices ar e compared to the single stellar population 
models of Schiavon (20071. These models include the 
effects of multiple element abundance ratios and allow es- 
timates of [Fe/H], [Mg/Fe], [C/Fe], [N/Fe], and [Ca/Fe], 
as well as mean lumino sity-weighted stellar age. In 
Graves & Schiavon ( |2008[ |, we created an IDL code pack- 
age, EZ_Ages, which automates the process of determin- 
ing stellar population mean age, [Fe/H], and abundance 
ratios for a set of Lick index measurements. We use 
EZ_Ages with the "standard set" of Lick indices (H/3, 
(Fe) 7 , Mg b, C 2 4668, CNi, and Ca4227) to determine 
ages, [Fe/H], and abundance ratios for the 84 stacked 



spectra in our sample with adequate S/N (see ^3.1) 8 . 
A detaile d description of the mode lling process is pre- 
sented in Graves fc Schiavon] ( |2008[ ) , as well as rigorous 
testing of the method on Galactic globular clusters and 
a co mparison with s tellar population modelling results 



from |Thomas et al. ( 2005 1 . The interested reader is re- 
ferred to that work for more information on the mod- 
elling process. It is worth noting here, however, that 
the models are computed at fixed [Fe /H] , rather than at 
fixed total metallicity ([Z/H]) as in some stellar popu- 
lation models. Total metallicity is dominated by oxy- 
gen, which is unmeasurable in unresolved stellar pop- 
ulations with current techniques. In this analysis, we 
specifically discuss [Fe/H] and [Mg/H], elements which 
can be measured from low-resolution spectra, and avoid 

7 (Fe) = | (Fe5270 + Fe5335) 

8 For stacked spectrum #76, the H/3 line was weaker than any of 
the models. We therefore computed abundances for this spectrum 
using H7i? instead of H/3. Bec ause Hyf produces syste matically 
younger age estimates than H/3 (IGraves & Schiavon|2008| , we have 
excluded this spectrum from the age analysis while retaining it for 
the abundance analysis. 
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Fig. 4. — Stellar population modelling results, showing mean luminosity-weighted stellar age, [Fe/H], [Mg/Fe], and [Mg/H] as a function 
of galaxy cr. Black (gray) points and error bars indicate high (low) S/N measurements and their associated median errors. Dashed lines 
show linear least squares fits of the stellar population properties as function of cr, based on the high S/N (black) data only. Mean stellar 
age, [Fe/H], [Mg/Fe], and [Mg/H] all increase with increasing cr. Age and [Fe/H] both show substantial scatter at fixed a, with total spread 
4—5 times the expected spread due to measurements errors, indicating genuine underlying population variations at fixed cr. [Mg/H] and 
[Mg/Fe] show less scatter, only ~ 2 times that expected due to measurement errors. The [Mg/H]— a relation is particularly strong and 
tight, nearly consistent with measurement errors, particularly at the high— cr end. 



making statements about [Z/H]. 

Figure [3] shows the measured values of (Fe) and H/3 
for t he stacked spectr a plotted over a set of model grids 
from Schiavon (20071. Data points are color-coded by 
a as indicated in the figure key. Crosses, filled circles, 
and open triangles represent the low-SB, midplane, and 
high- SB data, respectively. The various bins in R e are 
not indicated. Two trends emerge clearly from the data. 
As a increases, galaxies tend to have stronger (Fe) and 
weaker H/3. At fixed cr, as AI e increases, galaxies tend 
to have stronger (Fe) and stronger H/3. We will discuss 
the detailed implications of these trends in Q here they 
are merely presented to illustrate the linestrength data. 

The solid lines in Figure [3] show models of constant 
[Fe/H], from -0.7 to +0.2, as labelled. The dotted lines 
show models of constant age, from 2.2 to 14.1 Gyr. As 
can be seen from the orientation of the grid lines, H/3 
is predominantly sensitive to stellar population age such 
that younger populations have stronger H/3, while (Fe) is 
dominated by [Fe/H] such that more Fe-rich populations 
have stronger (Fe) . The effects of the age-metallicity de- 
generacy are apparent in the fact that H/3 is (mildly) sen- 
sitive to [Fe/H] and (Fe) is sensitive to age. Because of 
the degeneracy, there is no one-to-one mapping between 
H/3 and age, or (Fe) and [Fe/H]; the grid lines are not per- 
fectly vertical and horizontal. However, the model grids 
in these particular indices are close enough to orthogonal 
that they serve to break the age-metallicity degeneracy. 



Other comb inations of indices (e.g., H fcp and Fe4383, see 
Figure 9 of Graves & Schiavon 2008 ) are more affected 
by the degeneracy, causing the model grids to collapse 
down on top of one another and making them less ro- 
bust indicators of age and metallicity. 

EZ_Ages computes statistical errors in each of the stel- 
lar population parameters, determined from the mea- 
surement errors in the Lick indices. The black cross 
in the upper left corner of Figure [3] illustrates the me- 
dian errorbars in the line strenth measurements. The 
non-orthogonality of the model grids due to the age- 
metallicity degeneracy causes the derived errors to be 
correlated. For example, an overestimate of (Fe) will re- 
sult in both a higher derived value of [Fe /H] and a slightly 
lower derived value for the age. The age-metallicity de- 
generacy therefore causes correlated errors in the stel- 
lar population parameters, which must be taken into ac- 
count when interpreting results. In general, however, 
the very high S/N of the stacked spectra in this analy- 
sis make the stellar population modelling results robust 
to such effects. The effect o f cor related errors will be 
discussed in greater detail in j ]4.3| 

This paper discusses only the mean luminosity- 
weighted ages, [Fe/H], the abundance ratio [Mg/Fe], and 
[Mg/H] (which is simply [Mg/Fe] + [Fe/H]). The re- 
sults of the stellar population fits are given in Table 2, 
along with the Lick indices from which these parameters 
are determined. Discussion of other abundance ratios 
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Fig. 5. — Stellar population modelling results, showing mean luminosity-weighted stellar age, [Fe/H], [Mg/Fe], and [Mg/H] as a function 
of galaxy R e . Black (gray) points and error bars indicate high (low) S/N measurements and their associated median errors. No strong 
correlations are observed between R e and any stellar population properties, implying that galaxy star formation histories are nearly 
independent of galaxy size. 



is deferred to future work. An important caveat about 
the stellar population parameters discussed here is that 
they represent mean, luminosity-weighted averages over 
all the stars in a galaxy, which additionally have been 
averaged over all galaxies which contribute to a stacked 
spectrum. A younger mean age measured in one stacked 
spectrum compared to another does not make a defini- 
tive statement about the ages of all stars in all galaxies 
in that bin: it could also mean that all of the galax- 
ies in that bin contain sub-populations of younger stars, 
or that some of the galaxies in the bin have generally 
younger stars. There is also a strong degeneracy between 
the age of a younger subpopulation and the mass frac- 
tion it comprises: mean ages can be skewed younger by a 
small population of very young stars, or by a larger pop- 
ulation of somewhat less young stars. Age measurements 
should therefore be treated as a statistical description of 
galaxies in the bin, rather than as "true" ages for all the 
stars in an individual galaxy. 

4. STELLAR POPULATIONS IN FUNDAMENTAL PLANE 

SPACE 

4.1. Age, [Fe/H], [Mg/H], and [Mg/Fe] as Functions of 
a, R e , and I e Separately 

We begin by examining the variations in stellar popu- 
lation properties with each of the FP parameters a, R e , 
and I e individually. 

Figure [4] shows the EZ_Ages results for mean 
luminosity- weighted age, [Fe/H], [Mg/H], and [Mg/Fe] 
as a function of the median log a in each galaxy bin. We 
have divided the data into u low-S/N" and "high-S/TV" 



subsamples in order to indicate the range of error values, 
where low-S/N data points are those whose statistical 
error estimates are > 0.015 dex higher than the median 
error value for the data as a whole. The black data points 
and error bars indicate the values associated with high- 
S/N data. Low-S/N data are indicated in gray, along 
with the appropriate median error bar. Age, total Fe 
and Mg abundances, and the [Mg/Fe] abundance ratio 
all increase with increasing a. These relations are very 
similar to those demonstrated in Paper I, where galax- 
ies were binned first by a and then by L and color at 
fixed a. Dashed lines indicated linear least squares fits 
of the stellar population parameters as functions of a, 
computed using only the high S/N data (i.e., excluding 
the gray points). 

The age-cr and [Fe/H]-<r relations both show substan- 
tial scatter. Quantitatively, the spread around the mean 
relations (dashed lines) is ^4-5 times larger than the 
scatter expected due to the statistical errors in both age 
and [Fc/H], indicating that there is substantial spread 
in the star formation histories of galaxies at fixed a. 
The [Mg/Fe]-er and [Mg/H]-c relations show less scatter, 
with the observed spread representing only ^2 times the 
expected statistical scatter. The Mg abundance shows 
both the steepest correlation with a and the least intrin- 
sic scatter about the mean relation as compared to the 
observational errors, indicating that [Mg/H] varies only 
slightly at fixed a . 

Unlike the strong stellar population trends observed 
with tr, neither R e nor I e shows strong correlations with 
stellar population properties. Figure [5] shows the results 



8 



Graves et al. 



1.3 " 



u 1.0 - 



O 



0.8 - 



0.6 " 



2 



0.5 r 
0.4 r 
0.3 '- 
0.2 r 
0.1 r 
0.0 f- 
-0.1 r 



»* ** B 

- o" o o ft 

* ***** • 

* * 



■F 



i _ 



Ti 



2.0 2.5 3.0 3.5 
log / e (io/pc S ) 



— i — i — i — i — | — i — i — i — i — | — i— i— 

© 

***** 

* J&VvV 

8 * 



I : 



*© # *« <^ 
<? . ©o 0* f 

a -ft 







hi 



o.i 
o.o 

-0.1 „ 

X 

-0.2 > 
b. 

-0.3 '— , 

-0.4 

-0.5 

0.3 
0.2 

0.1 s 

0.0 w 
3 

-0.1" 

-0.2 

-0.3 



2.0 2.5 3,0 
log I. (L /pc e ) 



3.5 



Fig. 6. — Stellar population modelling results, showing mean luminosity-weighted stellar age, [Fe/H], [Mg/Fe], and [Mg/H] as a function 
of galaxy I e . Black (gray) points and error bars indicate high (low) S/N measurements and their associated median errors. Weak trends 
are visible, such that stellar population age decreases slightly with increasing I e , while [Fe/H] and [Mg/H] increase slightly with increasing 
I e , but none of these trends are compelling. [Mg/Fe] appears entirely independent of I e . 



of the stellar population modelling as a function of the 
median logi? e for each galaxy bin. There are no strong 
trends with R e in any of the properties modelled here. 
There is some evidence for slight increases in [Fc/H] and 
[Mg/H] with increasing R e , but the dependence is weak 
and the scatter is large. Also, at large R e there do not 
appear to be any "young" galaxies with mean ages < 7 
Gyr. However, it is evident from Figure [2] that there are 
very few galaxies with large R e and low a. Since Figure [4] 
showed that low-c galaxies are most likely to have young 
ages, the lack of "young" galaxies with large R e may be 
due to the way galaxies populate the FP, rather than any 
genuine relation between larger R e and older ages. 

Finally, Figure :G: shows age, [Fc/H], [Mg/Fe], and 
[Mg/H] plotted against the median log I e for each galaxy 
bin. Again, there are no strong correlations between 
log/ e and any of the stellar population parameters. 
There is some indication of weak correlations between 
log/e and [Fe/H] and [Mg/H] such that the high-SB 
galaxies have higher Fe and Mg abundances. There is 
also some evidence for a weak anti-correlation between 
log/ e and age such that low-SB galaxies have older ages, 
but none of these trends is compelling. The [Mg/Fe] 
abundance ratio is completely flat with log/ e . 

A clear result from Figures [4j|6] is that stellar popula- 
tions on the FP vary primarily with a. Neither R e nor I e 
appears to be closely related to stellar population prop- 
erties. However, the lack of general trends with R e and 
I e may hide weak correlations with residuals from the 
strong a relationships. Although Figure [4] shows clear 
correlations with a in all stellar population parameters, 



there is some scatter at fixed a, particularly in age and 
[Fe/H]. The next section examines how stellar popula- 
tions vary throughout the FP as a function of all three 
FP parameters simultaneously 

4.2. Mapping Stellar Populations Along and Through 
the Fundamental Plane 

By binning galaxies along and through the thickness of 
the FP, then stacking their spectra and performing stel- 
lar population modelling on the high S/N average spec- 
tra, we have measured the typical age, [Fe/H], [Mg/H], 
and [Mg/Fe] for galaxies at each point on, above, and 
below the FP. Each of these properties can be mapped 
throughout the 3-D space of the FP, thereby associating 
information about typical star formation histories with 
the structural properties of galaxies. This is a somewhat 
challenging exercise in data visualization, as we are try- 
ing to understand the variations and co- variations of four 
different stellar population parameters throughout a 3-D 
structure space. This effectively involves probing a 7- 
D parameter space, although by focusing on one stellar 
population parameter at a time, we can limit it to a 4-D 
space. In this and the next section, we attempt to codify 
trends in this 7-D space and understand how the various 
galaxy properties correlate with one another. We begin 
by mapping each of the four stellar population properties 
along and through the FP. 

FigurelTlshows the results for mean luminosity-weighed 
age. The 'figure is constructed to be similar to Figure [2] 
which illustrated the bin definitions for the 3-D FP space. 
The three lower panels show the three slices in A log I e , 
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Fig. 7. — Mean luminosity-weighted stellar population age, mapped across the FP. The three lower panels show the three slices in surface 
brightness, as in Figure|2] Within each slice, galaxy a and R e are plotted, with gray points sh owing individual galaxies in the sample. The 
overplotted color contours show the typical stellar population age at each point on the FP (see j ]4.2| for details). The dashed lines show lines 
of constant Md yn , assuming Af t j yn <x cr 2 R e . Within each slice, age increases with increasing <r, as seen in Figure[5] Lines of constant age 
run approximately vertically, indicating that mean stellar population age is independent of R e at fixed a. The upper panels show stellar 
population age as a function of a for each FP slice. Different values of R e are indicated by different line styles. Comparing the typical 
ages between the various slices in I e , it is clear that galaxies with lower surface brightness than the FP (left panels) have older ages than 
those on the FP midplane (center panels), while galaxies with higher surface brightness than the FP (right panels) have younger ages, at 
the same value of a and R e . Despite the lack of strong stellar population trends with absolute galaxy I e (c.f., Figure[6jl, stellar population 
age at fixed a does appear to depend on I e residuals from the FP, such that low-SB (high-SB) galaxies are older (younger) than those on 
the FP midplane. 



with the center panel showing the midplane slice through 
the FP, while the left and right panels show respectively 
the low-SB and high-SB slices from Figure [T] Within 
each FP slice, log a and log R e are plotted for the individ- 
ual galaxies as gray points (as in Figure |2| . Dashed lines 
show lines of constant Md yn . The color overlay indicates 
the mean stellar population age at each point in the pa- 
rameter space, with numerical values given by the color 
bar in the top of the leftmost panel. The age overlay 
was constructed as follows: for each stacked spectrum, 
the median values of log a and log R e are plotted for the 
galaxies in the bin, color-coded by the age measured in 
the stacked spectrum. Ages are then interpolated be- 
tween the median points to create a continuous map of 
stellar population age across the log er-log R e diagram. 
High-S/N and low- 5/ AT data are all included in the in- 
terpolation. 

The upper three panels of Figure [7J show mean 
luminosity-weighted age as a function of a for each FP 
slice. Different bins in R e are indicated by different line 
styles. 

Look first at the central panels of Figure [JJ which show 
the midplane of the FP. The trend from Figure [4] is vis- 
ible: galaxies with larger a are older (red) than galax- 
ies with lower a (blue). Although there are significant 



changes in age as a function of a, there seems to be lit- 
tle if any variation in age as a function of R e . Lines of 
constant age run roughly vertically in the diagram. This 
is consistent with Figure [5] which showed no clear sys- 
tematic variation in stellar population age as a function 
of i? e , independently of the other FP parameters. Here 
we see directly that, at fixed tr, there is no substantial 
dependence of stellar population age on R e . If the best 
predictor of stellar population age were Md yn rather than 
a, lines of constant age should follow the dashed line of 
constant Md yn - The near-verticality of the age contours 
indicates that a is a better predictor of stellar population 
age than is M dyn . 

If we look now at the galaxies from the low-SB (left 
panels) and high-SB (right panels) slices above and below 
the FP, we see essentially the same trends: age increases 
systematically with increasing er, but lines of constant 
age run nearly vertically, indicating that stellar popu- 
lation age is independent of R e at fixed a. However, 
comparing the age ranges (indicated by the color scale) 
between the different panels, there are systematic differ- 
ences. These are illustrated clearly in the upper panels. 
Galaxies in the low-SB slice of the FP (left panels) are 
systematically older than galaxies in the midplane of the 
FP (center panels) , which are in turn systematically older 
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Fig. 8. — Mean stellar population [Fe/H], mapped across the FP. The three lower panels show the three slices in surface brightness, as 
in Figure [5] Within each slice, galaxy a and R e are plotted, with gray points showing individual galaxies in the sample. The overplotted 
color contours show the typical [Fe/H] at each point on the FP. The dashed lines show lines of constant M l i yn , assuming Md yn oc o 2 R e . 
Within each slice, [Fe/H] increases with increasing <r, as in Figure H] [Fe/H] appears to be relatively independent of R B , with contours of 
equal [Fe/H] running nearly vertically in all slices. The upper panels show [Fe/H] as a function of a for each FP slice. Different values of 
R e are indicated by different line styles. Comparing the various slices in surface brightness, low-SB galaxies (left panels) appear to have 
lower [Fe/H] than galaxies on the midplane of the FP (center panels) at the same a and R e , while high-SB galaxies (right panels) have 
higher [Fe/H]. Combining these results with those apparent in Figure [7] suggests that, at fixed <r, stellar population age and [Fe/H] are 
anti- correlated. 



than galaxies in the high-SB slice (right panels). There 
are substantial overlaps in the age ranges in the three 
panels, but at fixed a the low-SB galaxies are older and 
the high-SB galaxies are younger than the typical galax- 
ies in the FP. This means that, although there are no 
clear trends in stellar population age with I e alone (Fig- 
ure [6]) , typical galaxy ages do vary with A log I e at a fixed 
point i n a an d R e . This is similar to the resu lts of |Forbes 



et al.| (|1998| and |Terlevich fc Forbesj ( |2002|) who found 
that surface brightness residuals from the FP correlate 
with stellar population age, although they explored only 
the general trend over all galaxies, not as a function of 
cr and R e . Our higher-resolution map is made possible 
through the large sample of 16,000 galaxies in this study. 

A similar map for [Fe /H] is shown in Figure [8j As 
seen in Figure El [Fe/H] increases systematically with 
a. Like stellar population age, [Fe/H] appears to de- 
pend weakly or negligibly on R e — contours of constant 
[Fc /H] are roughly vertical rather than following the lines 
of constant Md yn . Also similar to the age variations, 
there appear to be systematic variations in [Fe/H] be- 
tween the different slices in AlogJ e . This trend, how- 
ever, runs opposite to the one seen in age: at fixed cr, the 
low-SB galaxies (left panels) tend to be Fe-poor and the 
high-SB galaxies (right panels) tend to be Fe-rich. The 
general conclusion is that, not only do age and [Fe/H] in- 
crease with cr, they also vary systematically with A log I ei 



with low-SB galaxies showing older mean ages and lower 
[Fe/H], while high-SB galaxies show younger mean ages 
and higher [Fe/H]. This implies that, at fix ed cr, age and 
[Fe/H] are onij-correlated, as first shown by Trager et al. 
(20001. Section 4.3 will directly examine the correia- 



tions between various stellar population residuals from 
the mean trends with cr. 

Figure 9] shows the map for [Mg/H]. As expected from 
Figure [4 [Mg/H] varies strongly with cr. Like stellar 
population age and [Fe/H], [Mg/H] shows no depen- 
dence on R e . Unlike age and [Fe/H], [Mg/H] varies only 
mildly with Alog/ e - To the extent that [Mg/H] does 
vary at fixed a between the different panels, it follows 
the same trend as [Fe/H], with low-SB galaxies show- 
ing slightly lower [Mg/H] and high-SB galaxies showing 
slightly higher [Mg/H] than their counterparts in the 
midplane FP slice. However, the variation in [Mg/H] 
with Alog/g is only about half than that seen in [Fe/H]. 
FigureHshowed that [Mg/H] has the strongest and tight- 
est trend with cr. This 4-D map both confirms that there 
is little spread in [Mg/H] at fixed cr and illustrates that 
the small existing spread is correlated with Alog/ e . 

Finally, Figure 10 shows the map of [Mg/Fe] across 
the FP slices. As with the other stellar population pa- 
rameters, [Mg/Fe] depends strongly on cr and is indepen- 
dent of R e . Comparison between the panels shows that 
it also varies with Alog/ e . Unlike [Fe/H] and [Mg/H], 
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Fig. 9. — Mean stellar population [Mg/H], mapped across the FP. The three lower panels show the three slices in surface brightness, as in 
Figure[2] Within each slice, galaxy a and R e are plotted, with gray points showing individual galaxies in the sample. The overplotted color 
contours show the typical [Mg/H] at each point on the FP. The dashed lines show lines of constant Md yn , assuming M^ yn <x a R e . Within 
each slice, [Mg/H] increases strongly with increasing a, as in FigureH] Like age and [Fe/H], [Mg/H] appears to be relatively independent 
of R e , with contours of equal [Mg/H] running nearly vertically in all slices. The upper panels show [Mg/H] as a function of a for each FP 
slice. Different values of R e are indicated by different line styles. Comparing the various slices of the FP, there appears to be little variation 
in [Mg/H] between the low-SB galaxies (left panels), the FP midplane galaxies (center panels), and the high-SB galaxies (right panels) at 
fixed <t and R e . There is a slight tendency for high-SB galaxies to have higher [Mg/H], but the trend with Alogie is substantially weaker 
than that observed for [Fe/H] in Figure RS] 



both of which increase with increasing surface brightness, 
[Mg/Fe] decreases on average in high-SB galaxies (right 
panels) and increases on average in low-SB galaxies (left 
panels). In this respect, [Mg/Fe] behaves similarly to 
stellar population age, showing the highest values in low- 
SB, older galaxies. 

It is interesting that the low-SB slice of the FP 
(left panels) shows substantially more total variation in 
[Mg/Fe] than do the other panels. This trend is consis- 
tent through all three slices. Indeed, the high-/ e slice 
(right panels) shows only a small range of [Mg/Fe], while 
the low-ie slice (left panels) shows substantial variation. 
An equally valid way to state this (although perhaps 
harder to see clearly in Figure [To| would be to say that 
high-cr galaxies (at the right sideof each map) show sub- 
stantial variation in [Mg/Fe] between galaxies, but that 
the low-er galaxies (at the left side of each map) all have 
similarly low values of [Mg/Fe]. 

Taken together, these maps of stellar population pa- 
rameters along and through the FP show that age, 
[Fe/H], [Mg/H], and [Mg/Fe] all vary strongly with a, 
but that they also appear to vary with Alog/ e ([Mg/H] 
varies only mildly, while the other three parameters vary 
substantially). None of them varies significantly with 
R e . From these maps, it is clear that the spread around 
the a relations shown in Figure [4] is not just due to sta- 
tistical errors but in fact is systematic behavior that is 



correlated with surface brightness differences. This fur- 
thermore implies that there are substantial correlations 
between stellar population properties at fixed a. The 
next section looks explicitly at the correlations between 
residuals from the er-driven relations of Figure [4] in order 
to understand the systematic co- variation of the various 
stellar population properties. 

4.3. Residuals from the Stellar Population-cr Relations 

Stellar population residuals A (log age), A [Fe/H], 
A [Mg/H], and A [Mg/Fe] are defined with respect to the 
mean relations versus a by subtracting the linear fits 
shown in Figure H (dashed lines). Figure 11 shows var- 
ious correlations between these residuals. Filled circles 
show the high-S/iV data (corresponding to black points 
in Figure |4| , while empty circles show the low-S/N data 
(gray points in Figure [4]) . Data points are color-coded 
by Alog/ e , with blue representing the high-SB slice of 
the FP, green representing the midplane of the FP, and 
red showing the l ow-S B slice of the FP. 

As discussed in §3.2[ the age-metallicity degeneracy re- 
sults in correlated errors in the stellar population mod- 
elling process. Correlated errors arise because all ab- 
sorption indices are sensitive to both age and metallicity, 
resulting in model grid lines that are not orthogonal in 
index-index space (see Figurep I . The error ellipses in the 



lower left corners of Figure [T 



indicate the slope of the 
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Fig. 10. — Mean stellar population [Mg/Fe], mapped across the FP. The three lower panels show the three slices in surface brightness, as 
in Figure [2] Within each slice, galaxy a and R e are plotted, with gray points showing individual galaxies in the sample. The overplotted 
color contours show the typical [Mg/Fe] at each point on the FP. The dashed lines show lines of constant M^ym assuming M^yn ok cr 2 R e . 
Within each slice, [Mg/Fe] increases strongly with increasing cr, as in Figure [4] Like the other stellar population parameters, [Mg/Fe] 
appears to be relatively independent of R e , with contours of equal [Mg/Fe] running nearly vertically in all slices. The upper panels show 
[Mg/Fe] as a function of a for each FP slice. Different values of R e are indicated by different line styles. On average, low-SB galaxies (left 
panels) have higher [Mg/Fe] than those on the FP midplane (center panels), which in turn have higher [Mg/Fe] than the high-SB galaxies 
(right panels). The low-SB galaxies show a strong increase in [Mg/Fe] with cr, while the high-SB galaxies show very little overall variation 
in [Mg/Fe]. Another way to state this trend is that, at low a (leftmost galaxies in each panel), all galaxies have relatively low values of 
[Mg/Fe], while at high a (rightmost galaxies in each panel), there is a substantial range in [Mg/Fe]. 



correlated errors from the stellar population modelling, 
as determined b y Monte Carlo simulations in |Graves c£] 
Schiavon (2008 see Figure 3 of that work). The solid 
black and dashed gray error ellipses correspond to the 
median errors for high- and low-S/N data, respectively, 
based on the error estimates for the stellar population 
parameters determined by EZ_Ages. 

From Figure [IT] it is clear that there are significant 
correlations between the residuals from the mean rela- 
tions with a. These correlations are strongly associated 
with A log I e and thus represent real variations in star 
formation histories for galaxies with the same a. Al- 
though the correlated residuals in panels a-c are in the 
same direction as the correlated errors from the stellar 
population modelling process, the residuals trends can- 
not be explained merely by the correlated errors, for two 
reasons. In all cases, the black error ellipses are far too 
small to account for the observed spread in residuals, so 
the degeneracies in the modelling process cannot produce 
the observed spread in the data. Even more conclusively, 
the residuals depend strongly on Alog/ e . This cannot 
be explained by measurement errors in the absorption 
line strengths propagated through the correlated errors 
in the stellar population modelling because the sense of 
the measurement errors should be random, not depen- 
dent on A log I e . Thus the known degeneracies of the 



modelling process cannot account for the observed corre- 
lation of residuals. The appendix further demonstrates 
that correlated systematic errors due to multi-burst stel- 
lar population models also cannot be responsible for the 
observed trends. 

Figure [TTfc indicates that there is a strong anti- 
correlation between stellar population age and [Fe/H] at 
fixed cr, such that high- SB galaxies tend to be younger 
and more Fe-rich than their low-SB counterparts at the 
same cr. The total spread is at least a factor of two 
(0.3 dex) in both age and Fe abundance, with A [Fe/H] 
oc — 0.7A(log ag e). A similar anti-correlation was noted 



by_ Trager et al. ( 2000 1 and more recently by Smith et al.| 
1 )2008 bp , whose error estimates for individual galaxies 
were small enough to allow them to claim a genuine anti- 
correlation in the parameters. The clear dependenc e o f 
age and [Fe/H] residuals on A log I e shown in Figure Hi 
lays to rest any doubts that the anti-correlation could be 
due to errors alone, even though the slope of the anti- 
correlation is similar to the expected correlation of ob- 
servational errors. In c ontrast to the clear dependence 
on Alog/ e shown here, Trager et al. (20001 did not find 
any dependence of the age- [Fe / HJ anti-correlation on I e 
but in fact found that there was no improvement in sta- 
tistical fits of age or [Fe/H] when I e was included with 
cr as one of the fitting parameters. This is probably be- 
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Fig. 11. — Residuals from the mean stellar population trends with a, color-coded by surface brightness residuals from the FP (Alog/ e ). 
The A value for each stellar population parameter is determined by subtracting from the measured value the mean relation given by the 
dashed line in Figure [4] In each panel, solid (open) circles and black (gray) 1-<t error ellipses represent high (low) S/N data and the 
corresponding measurement errors. Error ellipses show the direction of the correlated errors in the age and abundance determinations. 
Data are color-coded by galaxy surface brightness, with galaxies in the low-SB, midplane, and high-SB FP slices shown in red, green, and 
blue, respectively. Stellar population residuals from the trends with a are clearly correlated with one another and with Alog/ e such that, 
at fixed <t, low-SB galaxies tend to be older, Fe-poor, Mg-poor, and enhanced in [Mg/Fe] compared to galaxies in the FP midplane, while 
high-SB galaxies are younger, Fe-rich, Mg-rich, and with lower [Mg/Fe]. 



cause they were fitting against absolute I e rather than 
A log I e at fixed a and R e (recall that Figure [6] shows 
only very weak dependence of stellar population param- 
eters on absolute log/ e )- There is a message here for 
multi-parameter fits in general: to detect a relation ver- 
sus a residual in a higher-dimensional space, it is neces- 
sary to fit against that residual explicitly; the variation 
may remain concealed when using the parent coordinate 
alone. 

A similar anti-correlation between residuals in age and 
Mg abundance is apparent in Figure |llb . Again, the 
anti-correlation is strongly dependent on Alog/ e . The 
strength of the anti-correlation is weaker than the age- 
[Fe/H] trend at fixed a (A [Mg/H] cx -0.4A(logage)). 
This is consistent with the fact that Figure [4] shows less 
spread in [Mg/H] at fixed a than in [Fe/H]. In addition 
to the total abundance of Mg at fixed a, there are also 
systematic variations in the abundance ratio [Mg/Fe]. 
Figure 11; shows that [Mg/Fe] is anti-correlated with 
[Fe/H] atnxed a, so that Fe-poor galaxies are more Mg- 
enhanced than their Fe-rich counterparts. As with age, 
[Fe/H], and [Mg/H], residuals from the log tr-[Mg/Fe] re- 
lation are correlated with A log I e . Finally, Figure \TT\ i 
shows that residuals in age are positively correlated with 
residuals in [Mg/Fe], such that older galaxies are more 
Mg-enhanced than younger galaxies at the same a. The 
statistical errors in age and [Mg/Fe] are uncorrelated 
(i.e., the error ellipses in panel d are horizontal) be- 



cause, although age errors in the stellar population mod- 
elling process affect individual abundance measurements, 
abundance ratio determinations are robust to modest age 
errors. 

In summary, there are correlations between the resid- 
uals from all of the various stellar population parameter 
trends with a. These depend on Alog/ e such that high- 
SB galaxies have younger ages, higher [Fe/H], higher 
[Mg/H], and lower [Mg/Fe] than low-SB galaxies with 
the same a and R e . These trends are similar to those re- 
ported in Paper I, where more luminous galaxies (those 
with higher AlogLy) were similarly younger, more Fe- 
rich, more Mg-rich, and with lower [Mg/Fe] with respect 
to less luminous galaxies at the same a. 

In this context, one can ask whether AlogLy (from 
Paper I) or Alog/ e (from this work) is the better dis- 
criminator of differing stellar populations at fixed a . A 
"better" discriminator can be defined as the binning pa- 
rameter which maximizes age and [Fe/H] differences at 
fixed a (i.e., does the best job of separating galaxies 
with differing star formation histories from one another). 
Figure [12| compares the discriminatory power of binning 
galaxies by A log I e (as in this analysis) versus binning 
by AlogLy (as in Paper I). In each analysis, galaxies 
were divided into "bright" , "medium" , and "faint" bins 
at fixed a, based either on Alog/ e or on AlogLy. The 
"medium" brightness bins are not included in Figure [12] 
because the "bright" and "faint" bins best illustrate the 
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Fig. 12. — Comparison of the discriminatory power of AlogLv 
and A log I e for distinguishing different galaxy stellar populations 
at fixed a. Top: Black (gray) histograms show the age and 
[Fe/H] distributions for "brighter" ("fainter") galaxies, as iden- 
tified by A log I e . Galaxies with differing stellar populations are 
well-separated by the A log/,, criterion. Bottom: Black (gray) his- 
tograms show the age and [Fe/H] distributions for brighter (fainter) 
galaxies, as indcntified by A log Ly for the same sample of galaxies 
(data from Paper I). The AlogLv criterion also does a reason- 
able job of separating galaxies with different stellar populations, 
but overall provides less clear separation than A log I e : the galax- 
ies in each bin are more mixed (in age and [Fe/H]) and the total 
spread between the peaks (in both age and [Fe/H]) is smaller. Thus 
A log I e is preferable for distinguishing galaxies with different star 
formation histories at fixed <r. 



total spread in stellar population properties. 

The top panels of Figure [12] show histograms of the 
residuals from the age-a and[Fe/H] -a relations based 
on the Alog/ e used in this analysis. Sorting galaxies by 
Alog/e does an excellent job of separating out galax- 
ies with differing star formation histories: the galaxies 
with the high-SB (black histograms) have the youngest 
mean ages and the highest [Fe/H], while those with the 
low-SB (gray histograms) have the oldest mean ages and 
lower [Fe/H]. The separation between the two popula- 
tions is very clear. If AlogLy is used to bin galaxies at 
fixed a (lower panels), the results are generally similar. 
However, the differences between high and low peaks are 
larger in both cases with Alog/ e than with Alogiy 
This is as expected since A log I e is a "pure" stellar pop- 
ulation measure whereas A log Ly also contains informa- 
tion from Alogi? e (at fixed a). As illustrated in Figures 
[7]jT0| A log R e is irrelevant to stellar population param- 
eters, and thus its presence in A log Ly serves only to 
dilute the information in that quantity. To conclude, if 
one is searching for the best way to identify the youngest 
(or most Fe-rich) galaxies at a given a, A log I e is a bet- 
ter discriminant than Alogiy. There may be cases 
where radii are not measurable or not available for in- 
dividual galaxies, particularly in high redshift surveys, 
in which case A log Ly still provides reasonable discrim- 
inatory power. 

The stellar population maps discussed in §4. 2 1 showed 
very little variation in stellar population properties as a 
function of R e at fixed a. A direct comparison of stel- 
lar population residuals and their dependence on R e is 
shown in Figure |13| Here, stellar population residuals 
are shown as a function of A log R e , where A log R e is de- 



fined as the difference between the median log R e for each 
bin and the median log R e for all galaxies in the same a 
range. The weak dependence of R e on a has therefore 
been removed, highlighting the effect of R e alone. It is 
clear that the stellar population residuals arc unrelated 
to galaxy size at fixed a. 

This is an interesting and perhaps surprising result. 
The total Md yn scales as a 2 R e . Thus, if the funda- 
mental galaxy property driving stellar population vari- 
ations were galaxy mass, we would expect residuals from 
a trends to correlate with R e as well as with a. If this 
were the case, the contours of age, [Fe/H], [Mg/H], and 
[Mg/Fe] in Figures [7- 10 would follow the dashed lines of 
Mdyn- Instead, lines of constant population properties 
run nearly vertically in those figures, and the residuals 
shown here in Figure [13] are completely independent of 
R e . This suggests that it is in fact a and not total mass 
that is fundamentally related to the galaxy star forma- 
tion history. The emerging generation of cosmological 
models, including predictions for stellar velocity disper- 
sions, may be able to shed light on the different evolu- 
tionary histories of galaxies with the same total Md yn 
but different values of a (and vice versa). 

Taken together, these results indicate that a is the 
primary structural parameter controlling stellar popu- 
lation variations, and therefore the primary predictor of 
a galaxy's past star formation history. Indeed, a ap- 
pears to be more fundamental than total mass, because 
residuals from the mean trends between a and the var- 
ious stellar population properties do not correlate with 
R e . At fixed tr, there are additional, genuine correlations 
between stellar population parameters which seem to de- 
pend strongly on A log I e such that high-SB galaxies are 
younger, more Fe-rich, and more Mg-rich, but with less 
enhanced [Mg/Fe] than typical galaxies at the same a. 
Conversely, low-SB galaxies at the same a are older, Fe- 
poor, Mg-poor, and show significantly enhanced [Mg/Fe]. 
These correlations provide strong constraints on cosmo- 
logical models of galaxy formation and are an important 
tool for understanding past star formation histories of 
early-type galaxies. The implications for galaxy star for- 
mation histories will be discussed in detail in Paper III. 

4.4. Can the Observed Variation in Stellar Populations 
Be Caused by Measurement Errors in the FP 
Parameters ? 

We have argued that, because the thickness of the FP 
in I e is strongly correlated with stellar population varia- 
tions, that the thickness of the plane must be real, rather 
than merely scatter due to observational errors. But be- 
fore accepting this, it is worth checking that the observed 
stellar population trends cannot be caused by observa- 
tional errors in the FP parameters. It is obvious that, 
if the thickness of the FP in the I e dimension is due to 
errors in J e itself, we should not see any trends with I e . 
What about thickness in I e that is due to errors in R e 
and/or a that "scatter" galaxies into a different I e bin? 

The slope of the FP indicates that higher a galaxies 
tend to have higher I e , thus galaxies with observed high 
I e should be galaxies that have "scattered" in from high 
a (i.e., galaxies for which a has been underestimated). 
However, the galaxies scattered in from higher a would 
have older ages and higher [Mg/Fe], which is opposite to 
the observed trend with A log I e . It is therefore not possi- 
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Fig. 13. — Residuals from the mean stellar population trends with a as a function of residuals in R e . The residual AlogiJ e is computed 
from the median log R e value for each bin by subtract ing the median log R e for all galaxies at the corresponding cr (i.e., removing the weak 
dependence of R e on <r) . As expected from Figures [7|)10| there is no correlation between any stellar population residual and A log R e . 



ble that errors in a could produce the observed trends in 
I e and stellar population properties. In the case of R e , I e 
is calculated as I e — Lj2-KR\ so overestimating R e will 
result in underestimating I e . However, we have shown 
that stellar populations do not depend on R e and thus 
errors in measuring R e cannot explain the observed stel- 
lar population variations as a function of I e . This leads 
to the conclusion that the thickness of the FP in I e is in- 
trinsic to the galaxy population, rather than an artifact 
of measurement error, and therefore that the observed 
stellar population variation with A log I e is likewise real. 

5. DISCUSSION 

We have seen that, while stellar population properties 
and therefore galaxy star formation histories vary sub- 
stantially with cr, they do not vary with R e . Alt hough 
this is not a new result (cf., |Trager et al. 2000), it is 
somewhat perplexing, as it implies that cr, and not total 
galaxy mass, is the better predictor of galaxy star for- 
mation history. In the literature, a is often used as a 
proxy for galaxy mass in stellar population analyses; au- 
thors routinely refer to galaxies with higher a as "more 
massive" , implying that the two parameters are inter- 
changeable. Various authors have st udied stellar popu- 
lation variations as function of cr (e.g.,|Trager et al.|2000| 



Bernardi et ~aT1|2003b[ 
2008b), as functions of M, 
and as functions of L 



Thomas et"ap2005| ISmith et a 



Gallazzi et al. 2005 



e.g. , [Kuntschner fc Davies||199 
Thomas et al 



Terlevich fc Forbes|[2002l . |Thomas et al.| ( |2005p in fact 
convert between a and M* using the mean relation be- 
tween the two and present results for stellar population 
parameters as functions of either. 



Comparisons between all these works are somewhat 
confusing and do not, at first glance, appear to give con- 
sistent results, particularly when addressing age varia- 
tions in galaxies. This is in part because, although cr, 
M*, and L are all highly correlated, they are not exactly 
the same. It appears that stellar populations vary not 
only with cr, but that at fixed cr, they vary systemati- 
cally with the residuals from the mean a-L relation (as 
shown in Paper I) , so that stellar population trends with 
a look significantly different from trends with L. (Trends 
as a function of are very similar to trends with L, as 
M*/L varies little among these galaxies.) We argue in 
Paper I that the trends with cr are more fundamental 
and that tracking stellar population variations with L is 
not a good proxy for studying them as a function of a. 
The work presented here suggests that Md yn is also not a 
good proxy for cr and confirms that cr is the best predic- 
tor of galaxy star formation histories. Models of galaxy 
formation must reproduce this dependence on a. 

Simulations of galaxy mergers offer some insight into 
why R e might be unrela ted to galaxy stellar populations. 



Robertson et al. (20061 present a series of dissipational 
galaxy merger simulations, along with a study of how var- 
ious properties of the mergers affect the FP parameters 
of the resulting merger remnant galaxies. They find that 
the final R e of the merger remnant depends strongly on 
the orientation of the initial galaxy orbits and the total 
angular momentum of the pre-merger system. Different 
initial orientations produce a factor of two variation in 
the resulting R e of the merger remnant for re-simulations 
of collisions between identical progenitor galaxies. The 
stellar populations of a galaxy should be independent 
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of galaxy orbit angular momentum and orientation be- 
fore a merger; therefore dissipational mergers will tend 
to decouple the remnant R e from any original relation 
between galaxy size and star formation history should 
one exist. 

Subsequent dissipationless merging of early-type galax- 
ies will also tend to further scatter galaxies in R e , while 
having only a small effect on a and presumably having no 
effect on the stellar population of the galaxy. Simulations 
of dissipationless mergers betw een equal-mass systems 
by Boylan-Kolchin et al. (20051 result in remnants with 
R e around 70% (0.23 dexj higher than the pre-merger 
systems, while the a of the remnant differs by < 20% 
(0.06 dcx) from the pre-merger system. Thus equal-mass 
dissipationless merging modifies R e without create ma- 
jor changes in a or the stellar populations of galaxies, 
thereby scrambling any pre-existing correlations. 

These merger simulations suggest that any initial cor- 
relation between galaxy size and star formation history 
will be erased through extensive merging. Gas-rich dissi- 
pational mergers will produce remnants whose radii de- 
pend strongly on the orbital orientations of the progen- 
itors, scrambling any initial size-dependence. Variable 
quantities of subsequent gas-poor dissipationless merg- 
ers, which change R e but have only minor effects on a 
and no effect on the stellar population, will further serve 
to erase any original stellar population dependence on 

Re- 

6. CONCLUSIONS 

We have identified a sample of ~16,000 early-type 
galaxies from the SDSS spectroscopic Main Galaxy Sam- 
ple to study star formation histories in FP space. Bins 
of similar galaxies are defined along the FP in log a and 
log R e , and through the thickness of the FP in A log I e ■ 
Spectra of the galaxies in each bin are coadded to pro- 
duce a set of high S/N mean spectra that span the FP. 
The mean spectr a are compared with the stellar popu- 
lation models of Schiavon ( 20O7J using the automated 
IDL code EZ_Ages (Graves <k Schiavon 2008), which de- 
termines the mean luminosity- weighted age, abundances 
[Fe/H], [Mg/H], and the abundance ratio [Mg/Fe] for 
each stacked spectrum. This allows us to study the vari- 
ations of stellar populations throughout FP space. 

The star formation histories of early-type galaxies are 
found to form a two-parameter family. Interestingly, 
these two parameters correlate with a and with A log I e 
but not with R e . In other words, the two-parameter fam- 
ily of star formation histories vary across one dimension 
of the FP and through the thickness of the FP, but do 
not seem to vary with the second dimension across the 
FP. A detailed summary of our results follows. 

1. Stellar population age, [Fe/H], [Mg/H], and 
[Mg/Fe] all increase with increasing galaxy a. This 
implies a close relationship between a and the star 
formation history of a galaxy. Age and [Fe/H] show 
substantial spreads at fixed a, while [Mg/Fe] and 
[Mg/H] vary less. The relation between a and 
[Mg/H] is particularly strong and tight, showing 
very little spread at high a and only modest spread 
at low a. 

2. Stellar population properties are independent of 
galaxy size, R e - This suggests that a is a better 



predictor of the star formation history of a galaxy 
than total mass since dynamical mass scales as 
Mdyn oc a 2 R e - This can be understood in the con- 
text of merger-driven galaxy formation, where R e 
is strongly affected by the initial orbital conditions 
of the merging galaxies and may be subsequently 
further altered by dissipationless merging (which 
does not affect the stellar population). 

At fixed a and R e , stellar populations vary sub- 
stantially with surface brightness residuals from the 
FP, parameterized by A log I e . This has a num- 
ber of implications. Firstly, the thickness of the 
FP is real rather than merely due to measurement 
error. Furthermore, the thickness of the FP rep- 
resents an age sequence, such that galaxies with 
higher surface brightness than the midplane of the 
FP have younger ages, while galaxies with lower 
surface brightness than the midplane have older 
ag es. This result is i n agreement with the wor k 



of Forbes et al. 
and Treu et al.l 



Tcrlcvich & Forbes (2002) 



In addition to age variations through the thick- 
ness of the FP, there are corresponding variations 
in all other stellar population properties studied 
here. The variations in age, [Fe/H], [Mg/H], and 
[Mg/Fe] at fixed a are correlated with one another 
such that the galaxies with higher surface bright- 
ness than the FP midplane are younger, signifi- 
cantly more Fe-rich, only slightly more Mg-rich, 
and have lower [Mg/Fe] than their counterparts 
on the FP. Similarly, galaxies with lower surface 
brightness than the FP midplane are older, less Fe- 
rich, slightly less Mg-rich, and have substantially 
enhanced [Mg/Fe] relative to their counterparts on 
the FP. The general age-metallicity anti-correlation 
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Although these 

co-variations are in the directions of the correlated 
errors in stellar population modelling, neither the 
statistical errors nor the known systematic errors 
in the modelling process can explain the variations, 
strongly suggesting that they are real effects. 

5. The variations with A log I e at fixed a and R e 
demonstrated in this analysis are consistent with 
the results of Paper I, which showed similar stel- 
lar population variations with A log Ly at fixed a. 
However, the FP residual A log I e does a somewhat 
better job of distinguishing young Fe-rich popula- 
tions from old Fe-poor populations at the same a 
and is thus preferable to A log Ly for studying vari- 
ations in the star formation histories of early-type 
galaxies. 

This work demonstrates that the stellar populations 
of early-type galaxy form a two-parameter family. This 
two-parameter family maps onto a cross-section through 
the FP. Paper III in this series will explicitly quantify 
this mapping, as well as translating the derived stellar 
population parameters into model star formation histo- 
ries. 
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APPENDIX 



FROSTING MODELS 



Section [4~3l showed that correlated statistical errors cannot account for the co- variation of stellar population properties 
with A log I e at fixed a. The you nger ages observed for high-SB galaxies are genuine variations, in agreement with 



the results of |Forbes et al. (1998). However, the stellar population models used in this work are single-burst models. 
Because young stars are luminous and contribute disproportionately to integrated galaxy light and Balmer absorption 
line strengths, a small population of young stars can skew the "mean" age to younger values than a true, mass-weighted 
mean age. The effect is largest when the age difference between sub-populations is large compared to the mean stellar 
age, so that the youngest age measurements in our data are likely to be somewhat too young. Determinations of age 
and [Fe/H] from absorption indices are anti-correlated, so a young component would not only skew the mean age to 
younger values but should also skew the measured [Fc/H] to higher values. Because [Mg/Fe] is an abundance ratio. 



it is relatively insensitive to the correlated errors in age and [Fe /H] (see Figure 3 of Graves & Schiavon ( 2008 ) . The 
abundance [Mg/H] is calculated as [Fe/H] + [Mg/Fe], so errors in [Fe/H] will translate directly into comparable errors 
in [Mg/H]. 



Figure 14 shows absorption line strength predictions for H/3 and (Fe) from stellar population models. The black 
grid shows a set of single-burst models from |Schiavon| (|2007[) with solar abundance ratios. Solid lines indicate lines of 
constant [Fe/H], while dotted lines indicate lines of constant age, as labelled in the figure. In addition to the grid of 
single burst models, two different frosting models are shown. The dark gray line shows the effect of adding a burst 
with age 1.2 Gyr to an underlying stellar population with an age of 14.1 Gyr. Both the old and young population 
components have [Fe/H]= —0.4. Diamonds indicate the fraction of galaxy light contributed by the young burst, 
starting at the bottom at 0% and increasing in increments of 5% moving upward along the line. This line shows the 
trajectory in H/3-(Fe) space produced by adding young bursts at the same [Fe/H] to an underlying old population. 
Though this model is computed for a burst, the direction traced by adding a range of younger stars would be the same 
or tilted even more strongly to the left. The light gray line with "X" symbols shows a similar frosting model, but here 
the young 1.2 Gyr burst is four times more Fe-rich than the underlying old population, with [Fe/H]= +0.2. 

The colored points show the data from our stacked galaxy spectra. Different colors represent the different bins in 
cr, with purple indicating the lowest-cr bin and red indicating the highest-cr bin. At each value of cr, the three circles 
indicate the three bins in Alog/ e , with the smallest circle showing the low-SB bin and the largest circle showing the 
high-SB bin. The bins have been averaged over R e . 

It is interesting that the slopes between the low-SB and midplane bins (small and medium circles) are substantially 
different from those between the midplane and high-SB bins (medium and large circles). Comparing first the low-SB 
and midplane bins, it is clear that the differences between low- SB and midplane galaxies cannot be explained at all 
by a small frosting of younger stars at the same [Fe/H]. The same-Fe frosting model (dark gray line with diamonds) 
is nearly orthogonal to the change in H/3-(Fe) space between the low-SB and midplane bins. Even adding a frosting 
population with four times larger Fe abundance (light gray lines with exes) does not come close to explaining the (Fe) 
line strength differences between the low-SB and midplane galaxies. The low-SB galaxies must have genuinely lower 
values of [Fe/H] than the galaxies at the FP midplane, as well as having older ages. 

Line strengths differences between the midplane galaxies and the high-SB galaxies are more similar to the H/3-(Fe) 
trajectories expected for frosting models, particularly for low-cr galaxies. At high cr, neither the constant [Fe/H] nor 
the Fe-rich frosting models can account for the observed increases in (Fe) in the high-SB galaxies. For the lowest a 
bin, the difference between the midplane galaxies and the high-SB galaxies is consistent with a ~ 5% frosting of young 
stars at similar [Fe/H] to the underlying population — no genuine difference in [Fe/H] is needed to explain the data. 
However, the difference in [Fe/H] inferred from assuming a single burst population is only 0.04 dex. This indicates 
that the simplifying assumption of a single burst population only biases the estimated [Fe/H] by a small amount for 
the high-SB, lowest-cr galaxies, where the bias is strongest. This small bias in the measured [Fc/H] will translate into 
a similarly small bias in [Mg/H] . 

In general, a small frosting of younger stars, while adequate to explain the observed differences in stellar population 
ages at fixed cr, cannot also produce the observed variation in [Fe/H] and [Mg/H]. The variation in (Fe) linestrength 
at fixed cr requires a real variation in [Fe/H]. Thus the variations in observed abundance patterns at fixed a must be 
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Fig. 14. — Comparing observed index variations with two-burst "frosting" stellar population models. The black grid shows solar-ratio 
single burst models from [Schiavon (2007). Solid lines indicate lines of constant [Fe/H], while dotted lines indicate lines of constant age, as 
labelled. The dark gray trajectory shows the effect of adding a 1.2 Gyr burst with [Fe/H] = -0.4 to an underlying 14.1 Gyr old population 
with the same [Fe/H]. Diamonds indicate the fraction of the galaxy light at ~ 5000 A contributed by the young population, starting from 
0% at the bottom and increasing in increments of 5% to 45% at the top of the diagram. The light gray trajectory shows the effect of adding 
a 1.2 Gyr burst with four times larger Fe abundance ([Fe/H] = +0.2) to the underlying 14.1 Gyr population. "X"'s indicate the fraction 
of light contributed by the young Fe-rich population in increments of 5%. The indices measured in our stacked SDSS spectra are shown as 
colored points, with purple, blue, green, gold, orange, and red indicating the lowest through highest a bins, respectively. Small, medium, 
and large circles indicate the low-SB, midplane, and high-SB bins in Alog/ e , respectively. Differences between the low-SB and midplane 
bins cannot be produced merely by adding a young subpopulation of stars — a genuinely higher [Fe/H] population is required. Differences 
between the midplane and high-SB bins are compatible with a young subpopulation at the same [Fe/H] for the lowest a bin only — higher 
a bins also require a genuinely higher [Fe/H] population. The observed anti-correlation between mean age and [Fe/H] at fixed a therefore 
reflects genuine differences in the stellar population abundances and cannot be caused solely by systematic effects from using single burst 
stellar population models. 



due to genuine abundance differences between the galaxies, rather than an artifact of the modelling process for the 
younger galaxies. 
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Properties of Galaxy Bins and Stacked Spectra 
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0.42 


0.15 


2.89 


106 


349 


70 




0.5 < Re < 0.7 


-0.3 < A7 e < 


-0.1 


2.32 


0.56 


-0.13 


2.42 


54 


186 


71 






-0.1 < Ale < 


0.1 


2.31 


0.59 


0.01 


2.53 


504 


675 


72 






0.1 < A/ e < 


0.3 


2.31 


0.61 


0.14 


2.65 


161 


422 


73 




0.7 < Re< 0.9 


-0.3 < A/ e < 


-0.1 


2.33 


0.75 


-0.12 


2.19 


11 


89 
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-0.1 


< 


AI e 


< 


0.1 


2.32 


0.77 


0.03 


2.33 


213 


442 


75 




0.1 


< 


AI e 


< 


0.3 


2.32 


0.77 


0.14 


2.45 


67 


278 


76 


2.36 < a < 2.50 -0.1 < R e < 0.1 


-0.3 


< 


Ale 


< 


-0.1 


2.38 


0.02 


-0.15 


3.17 


19 


100 


77 




-0.1 


< 


Ale 


< 


0.1 


2.40 


0.03 


-0.03 


3.28 


37 


168 


78 




0.1 


< 


AI e 


< 


0.3 










3 


53 


79 


0.1 < R e < 0.3 


-0.3 


< 


Ale 


< 


-0.1 


2.39 


0.19 


-0.15 


2.94 


32 


149 


80 




-0.1 


< 


Ale 


< 


0.1 


2.38 


0.22 


-0.03 


3.03 


71 


266 


81 




0.1 


< 


AI e 


< 


0.3 










3 


58 


82 


0.3 < R e < 0.5 


-0.3 


< 


Ale 


< 


-0.1 


2.38 


0.40 


-0.13 


2.68 


35 


170 


83 




-0.1 


< 


Ale 


< 


0.1 


2.39 


0.42 


-0.01 


2.82 


149 


408 


84 




0.1 


< 


Ale 


< 


0.3 


2.37 


0.41 


0.14 


2.97 


12 


116 


85 


0.5 < Re < 0.7 


-0.3 


< 


Ale 


< 


-0.1 


2.38 


0.56 


-0.13 


2.49 


26 


152 


86 




-0.1 


< 


Ale 


< 


0.1 


2.39 


0.61 


0.01 


2.61 


242 


527 


87 




0.1 


< 


Ah 


< 


0.3 


2.39 


0.63 


0.13 


2.70 


43 


244 


88 


0.7 < Re < 0.9 


-0.3 


< 


Ah 


< 


-0.1 










3 


50 


89 




-0.1 


< 


Ale 


< 


0.1 


2.39 


0.76 


0.02 


2.42 


153 


420 


90 




0.1 


< 


Ah 


< 


0.3 


2.40 


0.77 


0.12 


2.54 


26 


190 



Note. — 7 e is computed in the V-band. 
* Quoted values arc in log units, i.e., logo - , logi? e , and A log I e . 
f Total number of galaxies in the stacked spectrum. 
| Effective median S/N of the stacked spectrum. 
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TABLE 2 

Lick Indices and Stellar Population Properties of Stacked Spectra 



H/3 


(Fe> 


Mg b Age 1 " 


[Fc/H] 


[Mg/H] [Mg/Fc] 


(A) 


(A) 


(A) (Gyr) 


(dex) 


(dcx) (dcx) 



1 1.82 ± 0.08 2.11 ± 0.09 3.29 ± 0.08 10.3 ± 1.5 -0.48 ± 0.05 -0.22 ± 0.08 0.26 ± 0.06 

2 2.00 ± 0.04 2.24 ± 0.05 3.10 ± 0.04 7.1 ± 0.5 -0.33 ± 0.03 -0.18 ± 0.03 0.15 ± 0.02 

3 2.23 ± 0.06 2.50 ± 0.07 2.93 ± 0.06 4.1 ± 0.5 -0.09 ± 0.04 -0.06 ± 0.05 0.03 ± 0.03 

4 1.81 ± 0.05 2.25 ± 0.05 3.21 ± 0.05 10.1 ± 0.9 -0.39 ± 0.03 -0.23 ± 0.03 0.16 ± 0.02 

5 1.99 ± 0.03 2.40 ± 0.03 3.27 ± 0.03 6.9 ± 0.3 -0.22 ± 0.01 -0.10 ± 0.03 0.12 ± 0.02 

6 2.20 ± 0.03 2.35 ± 0.04 3.02 ± 0.03 4.7 ± 0.3 -0.19 ± 0.02 -0.09 ± 0.03 0.10 ± 0.02 

7 1.86 ± 0.05 2.36 ± 0.05 3.30 ± 0.04 8.9 ± 0.7 -0.29 ± 0.03 -0.15 ± 0.04 0.14 ± 0.03 

8 2.07 ± 0.03 2.36 ± 0.03 3.27 ± 0.03 5.8 ± 0.3 -0.23 ± 0.02 -0.07 ± 0.03 0.16 ± 0.02 

9 2.17 ± 0.05 2.37 ± 0.05 3.37 ± 0.04 4.9 ± 0.4 -0.17 ± 0.03 0.02 ± 0.04 0.18 ± 0.03 

10 1.86 ± 0.09 2.48 ± 0.11 3.36 ± 0.09 8.8 ± 1.4 -0.22 ± 0.05 -0.14 ± 0.07 0.08 ± 0.05 

11 1.96 ± 0.05 2.45 ± 0.06 3.55 ± 0.05 7.2 ± 0.6 -0.19 ± 0.03 -0.03 ± 0.04 0.16 ± 0.02 

12 2.12 ± 0.08 2.27 ± 0.09 3.13 ± 0.08 5.4 ± 0.7 -0.26 ± 0.05 -0.09 ± 0.07 0.17 ± 0.05 

13 ... ± ± ± ± ± ± ± . . . 

14 1.78 ± 0.09 2.51 ± 0.11 3.40 ± 0.09 10.1 ± 1.6 -0.24 ± 0.05 -0.15 ± 0.07 0.09 ± 0.05 

15 ... ± ± ± ± ± ± ± . . . 

16 1.77 ± 0.06 2.49 ± 0.06 3.35 ± 0.06 10.4 ± 1.0 -0.26 ± 0.03 -0.18 ± 0.04 0.08 ± 0.03 

17 1.92 ± 0.03 2.42 ± 0.03 3.31 ± 0.03 7.8 ± 0.4 -0.23 ± 0.02 -0.12 ± 0.03 0.12 ± 0.02 

18 2.21 ± 0.04 2.46 ± 0.05 3.19 ± 0.04 4.5 ± 0.6 -0.11 ± 0.03 -0.02 ± 0.03 0.09 ± 0.02 

19 1.83 ± 0.04 2.42 ± 0.04 3.56 ± 0.03 9.2 ± 0.6 -0.25 ± 0.02 -0.08 ± 0.03 0.17 ± 0.02 

20 1.97 ± 0.02 2.47 ± 0.02 3.53 ± 0.02 6.9 ± 0.2 -0.16 ± 0.01 -0.01 ± 0.01 0.15 ± 0.01 

21 2.13 ± 0.03 2.43 ± 0.04 3.42 ± 0.03 5.2 ± 0.3 -0.15 ± 0.02 0.01 ± 0.03 0.16 ± 0.02 

22 1.82 ± 0.03 2.44 ± 0.04 3.41 ± 0.03 9.3 ± 0.6 -0.26 ± 0.02 -0.13 ± 0.03 0.13 ± 0.02 

23 1.93 ± 0.02 2.49 ± 0.02 3.57 ± 0.02 7.4 ± 0.2 -0.17 ± 0.01 -0.02 ± 0.01 0.15 ± 0.01 

24 2.10 ± 0.03 2.55 ± 0.04 3.49 ± 0.03 5.3 ± 0.3 -0.08 ± 0.02 0.05 ± 0.03 0.13 ± 0.02 

25 1.83 ± 0.06 2.32 ± 0.07 3.44 ± 0.06 9.4 ± 1.1 -0.32 ± 0.04 -0.13 ± 0.06 0.19 ± 0.05 

26 1.85 ± 0.03 2.51 ± 0.04 3.59 ± 0.03 8.6 ± 0.5 -0.19 ± 0.02 -0.05 ± 0.03 0.14 ± 0.02 

27 2.14 ± 0.06 2.50 ± 0.07 3.53 ± 0.06 4.9 ± 0.6 -0.09 ± 0.04 0.07 ± 0.05 0.16 ± 0.03 

28 1.96 ± 0.12 2.07 ± 0.15 3.28 ± 0.12 7.8 ± 1.6 -0.45 ± 0.09 -0.13 ± 0.12 0.32 ± 0.09 

29 1.85 ± 0.09 2.57 ± 0.10 3.57 ± 0.08 8.6 ± 1.3 -0.16 ± 0.05 -0.04 ± 0.07 0.12 ± 0.05 

30 ... ± ± ± ± ± ± ± . . . 

31 1.73 ± 0.04 2.46 ± 0.04 3.66 ± 0.04 11.1 ± 0.8 -0.27 ± 0.02 -0.10 ± 0.03 0.17 ± 0.02 

32 1.84 ± 0.02 2.49 ± 0.02 3.66 ± 0.02 8.7 ± 0.3 -0.20 ± 0.01 -0.04 ± 0.02 0.16 ± 0.02 

33 2.04 ± 0.03 2.51 ± 0.03 3.51 ± 0.03 6.0 ± 0.4 -0.13 ± 0.02 0.01 ± 0.03 0.14 ± 0.02 

34 1.80 ± 0.03 2.45 ± 0.03 3.79 ± 0.03 9.5 ± 0.5 -0.23 ± 0.02 -0.02 ± 0.03 0.21 ± 0.02 

35 1.86 ± 0.01 2.56 ± 0.01 3.79 ± 0.01 8.3 ± 0.2 -0.14 ± 0.01 0.03 ± 0.01 0.17 ± 0.01 

36 2.06 ± 0.02 2.61 ± 0.02 3.58 ± 0.02 5.5 ± 0.2 -0.05 ± 0.01 0.07 ± 0.02 0.12 ± 0.01 

37 1.74 ± 0.03 2.47 ± 0.03 3.85 ± 0.03 10.6 ± 0.5 -0.24 ± 0.02 -0.03 ± 0.03 0.21 ± 0.02 

38 1.84 ± 0.01 2.57 ± 0.01 3.86 ± 0.01 8.5 ± 0.2 -0.13 ± 0.01 0.05 ± 0.01 0.18 ± 0.01 

39 2.04 ± 0.02 2.60 ± 0.02 3.60 ± 0.02 5.8 ± 0.3 -0.06 ± 0.01 0.06 ± 0.02 0.13 ± 0.01 

40 1.76 ± 0.04 2.37 ± 0.05 3.67 ± 0.04 10.6 ± 0.8 -0.31 ± 0.03 -0.08 ± 0.04 0.23 ± 0.03 

41 1.81 ± 0.02 2.58 ± 0.02 3.84 ± 0.02 9.1 ± 0.3 -0.15 ± 0.01 0.02 ± 0.01 0.17 ± 0.01 

42 1.99 ± 0.03 2.65 ± 0.03 3.68 ± 0.03 6.4 ± 0.5 -0.04 ± 0.02 0.07 ± 0.03 0.11 ± 0.02 

43 1.89 ± 0.11 2.23 ± 0.13 3.59 ± 0.11 8.5 ± 1.6 -0.34 ± 0.07 -0.06 ± 0.09 0.28 ± 0.06 

44 1.72 ± 0.04 2.59 ± 0.05 3.89 ± 0.04 10.7 ± 0.8 -0.18 ± 0.02 -0.01 ± 0.03 0.17 ± 0.02 

45 1.84 ± 0.05 2.67 ± 0.06 3.91 ± 0.05 8.3 ± 0.8 -0.07 ± 0.03 0.08 ± 0.04 0.15 ± 0.03 

46 1.56 ± 0.04 2.46 ± 0.05 4.00 ± 0.04 15.2 ± 1.7 -0.19 ± 0.04 0.07 ± 0.04 0.26 ± 0.02 

47 1.75 ± 0.02 2.46 ± 0.02 3.96 ± 0.02 10.2 ± 0.4 -0.23 ± 0.01 0.02 ± 0.02 0.25 ± 0.02 

48 1.92 ± 0.03 2.61 ± 0.04 3.90 ± 0.03 7.3 ± 0.5 -0.08 ± 0.02 0.10 ± 0.03 0.18 ± 0.02 

49 1.69 ± 0.03 2.48 ± 0.03 3.99 ± 0.02 11.4 ± 0.5 -0.25 ± 0.01 0.01 ± 0.03 0.26 ± 0.03 

50 1.78 ± 0.01 2.61 ± 0.01 4.08 ± 0.01 9.4 ± 0.2 -0.12 ± 0.01 0.08 ± 0.01 0.21 ± 0.01 

51 1.97 ± 0.02 2.67 ± 0.02 3.85 ± 0.02 6.6 ± 0.2 -0.02 ± 0.01 0.12 ± 0.01 0.14 ± 0.01 

52 1.75 ± 0.03 2.50 ± 0.03 4.13 ± 0.02 10.2 ± 0.5 -0.20 ± 0.01 0.07 ± 0.03 0.27 ± 0.02 

53 1.76 ± 0.01 2.65 ± 0.01 4.10 ± 0.01 9.6 ± 0.2 -0.11 ± 0.01 0.08 ± 0.02 0.19 ± 0.02 

54 1.94 ± 0.02 2.67 ± 0.02 3.91 ± 0.01 7.0 ± 0.2 -0.03 ± 0.01 0.13 ± 0.02 0.16 ± 0.02 

55 1.75 ± 0.04 2.42 ± 0.04 4.20 ± 0.04 10.4 ± 0.7 -0.24 ± 0.02 0.08 ± 0.03 0.32 ± 0.02 

56 1.77 ± 0.01 2.66 ± 0.01 4.06 ± 0.01 9.4 ± 0.2 -0.10 ± 0.01 0.08 ± 0.01 0.18 ± 0.01 

57 1.90 ± 0.02 2.63 ± 0.02 3.91 ± 0.02 7.5 ± 0.3 -0.07 ± 0.01 0.10 ± 0.02 0.17 ± 0.01 

58 1.77 ± 0.07 2.34 ± 0.08 4.15 ± 0.07 10.1 ± 1.3 -0.29 ± 0.04 0.07 ± 0.06 0.36 ± 0.05 

59 1.70 ± 0.02 2.62 ± 0.03 4.05 ± 0.02 10.9 ± 0.5 -0.16 ± 0.01 0.03 ± 0.02 0.19 ± 0.02 

60 1.83 ± 0.04 2.57 ± 0.04 4.05 ± 0.04 8.5 ± 0.6 -0.12 ± 0.02 0.11 ± 0.03 0.23 ± 0.02 

61 1.59 ± 0.05 2.44 ± 0.06 4.35 ± 0.05 13.8 ± 1.3 -0.28 ± 0.03 0.05 ± 0.05 0.33 ± 0.04 

62 1.65 ± 0.03 2.59 ± 0.03 4.40 ± 0.03 11.8 ± 0.6 -0.17 ± 0.02 0.11 ± 0.03 0.28 ± 0.02 

63 1.76 ± 0.07 2.73 ± 0.09 4.25 ± 0.07 9.3 ± 1.1 -0.05 ± 0.04 0.14 ± 0.06 0.19 ± 0.05 

64 1.60 ± 0.03 2.59 ± 0.04 4.32 ± 0.03 12.9 ± 0.7 -0.19 ± 0.02 0.07 ± 0.03 0.26 ± 0.02 

65 1.66 ± 0.02 2.61 ± 0.02 4.40 ± 0.02 11.4 ± 0.3 -0.15 ± 0.01 0.13 ± 0.01 0.28 ± 0.01 

66 1.82 ± 0.03 2.79 ± 0.04 4.19 ± 0.03 8.4 ± 0.5 0.01 ± 0.01 0.17 ± 0.03 0.16 ± 0.02 

67 1.58 ± 0.03 2.59 ± 0.03 4.44 ± 0.03 13.3 ± 0.7 -0.19 ± 0.02 0.09 ± 0.03 0.28 ± 0.02 

68 1.70 ± 0.01 2.69 ± 0.01 4.34 ± 0.01 10.5 ± 0.2 -0.09 ± 0.01 0.14 ± 0.01 0.23 ± 0.01 

69 1.88 ± 0.02 2.73 ± 0.02 4.13 ± 0.02 7.5 ± 0.3 -0.00 ± 0.01 0.18 ± 0.02 0.18 ± 0.02 

70 1.60 ± 0.04 2.58 ± 0.04 4.47 ± 0.04 12.9 ± 0.9 -0.19 ± 0.02 0.11 ± 0.03 0.30 ± 0.02 

71 1.69 ± 0.01 2.69 ± 0.01 4.39 ± 0.01 10.7 ± 0.2 -0.10 ± 0.01 0.14 ± 0.02 0.24 ± 0.02 

72 1.87 ± 0.02 2.76 ± 0.02 4.20 ± 0.02 7.6 ± 0.2 0.01 ± 0.01 0.19 ± 0.01 0.18 ± 0.01 

73 1.60 ± 0.08 2.46 ± 0.09 4.68 ± 0.07 13.1 ± 1.8 -0.24 ± 0.05 0.17 ± 0.08 0.41 ± 0.06 

74 1.65 ± 0.02 2.67 ± 0.02 4.35 ± 0.02 11.5 ± 0.3 -0.12 ± 0.01 0.11 ± 0.02 0.24 ± 0.02 



23 



TABLE 2 

Lick Indices and Stellar Population Properties of Stacked Spectra 



75 


1 


.81 


± 





.02 


2 


.74 


± 





.03 


.3. 


22 


± 





.02 


8. 


.2 


± 





.4 


-0 


.05 


± 


0.01 


0. 


0.3 


± 


0. 


.82 


0. 


08 


± 





.82 


76 


1 


. 17 


± 


o 


.08 


2 


.50 


± 





.10 


1. 


59 


± 





.08 






± 




.t 


-0 


.19 


± 


0.06 


0. 


20 


± 


0. 


.08 


0. 


38 


± 





.05 


77 


1 


.62 


± 


o 


.04 


2 


.71 


± 


o 


.05 


4. 


58 


± 


o 


.04 


11 


.8 


± 


o 


.8 


-0 


.09 


± 


0.03 


0. 


18 


± 


o 


.04 


0. 


27 


± 


o 


.03 


78 






± 










± 










± 










± 










± 








± 










± 






79 


1 


53 


± 


o 


.05 


2 


.68 


± 





.06 


1. 


72 


± 





.05 


1 1 


.0 


± 


1 


.4 


-0. 


.11 


± 


0.03 


0. 


16 


± 


0. 


.05 


0. 


30 


± 





.04 


80 


1 


.62 


± 


o 


.03 


2 


.75 


± 





.03 


1. 


60 


± 





.03 


11 


.7 


± 


o 


.5 


-0. 


.07 


± 


0.02 


0. 


19 


± 


0. 


.03 


0. 


26 


± 





.03 


81 






± 










± 










± 










± 










± 








± 










± 






82 


1 


.55 


± 





.04 


2 


.55 


± 





.05 


1. 


67 


± 





.04 


13 


.7 


± 


1 


.1 


-0. 


.20 


± 


0.03 


0. 


17 


± 


0. 


.04 


0. 


37 


± 





04 


83 


1 


.58 


± 





.02 


2 


.71 


± 





.02 


1. 


67 


± 





.02 


12 


.1 


± 





.3 


-0 


.09 


± 


0.01 


0. 


19 


± 


0. 


.01 


0. 


28 


± 





.01 


84 


1 


.79 


± 





.07 


2 


.87 


± 





.08 


1. 


60 


± 





.06 


8. 


.4 


± 


1 


.0 


0. 


.05 


± 


0.03 


0. 


27 


± 





.04 


0. 


22 


± 





.03 


85 


1 


.55 


± 





.05 


2 


.66 


± 





.05 


4. 


63 


± 





.05 


13 


.5 


± 


1 


.2 


-0. 


.15 


± 


0.03 


0. 


15 


± 





.04 


0. 


30 


± 





.02 


86 


1 


.63 


± 





.01 


2 


.76 


± 





.02 


4. 


57 


± 





.01 


11 


.5 


± 





.3 


-0. 


.06 


± 


0.01 


0. 


19 


± 


0. 


.01 


0. 


25 


± 





.01 


87 


1 


.71 


± 





.03 


2 


.80 


± 





.03 


4. 


47 


± 





.03 


9 


.9 


± 





.5 


-0. 


.02 


± 


0.02 


0. 


19 


± 





.03 


0. 


21 


± 





.02 


88 






± 










± 










± 










± 










± 








± 










± 






89 


1 


.60 


± 





.02 


2 


.80 


± 





.02 


4. 


61 


± 





.02 


11. 


.9 


± 





.3 


-0. 


.05 


± 


0.01 


0. 


20 


± 


0. 


.01 


0. 


25 


± 





.01 


90 


1 


.74 


± 





.04 


2 


.83 


± 





.04 


4. 


36 


± 





.03 


9 


.1 


± 





.6 


0. 


.01 


± 


0.01 


0. 


19 


± 


0. 


.03 


0. 


18 


± 





.02 



Note. — f Ages arc determined from H/3 for all stacked spectra with the exception of spectrum 76, whose weak HP absorption falls outside 
the parameter space covered by the models. For this index, is used in the fitting process to determine abundances. For consistency, the age 

measure d for this spectrum is no t used in the analysis here, as ages derived from H-yp are typi can Y somewhat younger than those derived from 
Uf3 (see |Graves & Schiavon|2008[ . 



